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Abstract 



We show that there exist six parallel textures of the charged lepton and 
neutrino mass matrices with six vanishing entries, whose phenomenological 
consequences are exactly the same. These isomeric lepton mass matrices are 
compatible with current experimental data at the Su level. If the seesaw 
mechanism and the Fukugita-Tanimoto-Yanagida hypothesis are taken into 
account, it will be possible to fit the experimental data at or below the 2ct level. 
In particular, the maximal atmospheric neutrino mixing can be reconciled 
with a strong neutrino mass hierarchy in the seesaw case. 
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The recent solar [1], atmospheric [2], KamLAND [3] and K2K [4] neutrino oscillation 
experiments have provided us with very convincing evidence that neutrinos are massive 
and lepton flavors are mixed. In particular, the admixture of three lepton flavors involves 
two large angles 9i2 ~ 33° and 623 ~ 45° [5]. To interpret the observed bi-large lepton 
flavor mixing pattern, many phenomenological ansatze of lepton mass matrices have been 
proposed in the literature [6]. A very interesting category of the ansatze focus on texture 
zeros of charged lepton and neutrino mass matrices in a speciflc flavor basis, from which 
some nontrivial and testable relations between flavor mixing angles and lepton mass ratios 
can be derived. A typical example is the Pritzsch ansatz [7] of lepton mass matrices, 



in which six texture zeros are included and all non- vanishing entries are simply symbolized 
by x's. It has been shown in Ref. [8] that this ansatz can naturally predict a normal but 
weak neutrino mass hierarchy and a bi-large lepton flavor mixing pattern. If the seesaw 
mechanism is incorporated in the Fritzsch texture of charged lepton and Dirac neutrino 
mass matrices [9], one may obtain a similar flavor mixing pattern together with a much 
stronger neutrino mass hierarchy. 

The simplicity and predictability of Mi and Mj, in Eq. (1) motivate us to examine 
other possible six-zero textures of lepton mass matrices and their various phenomenological 
consequences. We flnd that there totally exist six parallel patterns of Mi and M^, with six 
texture zeros, as hsted in Table 1, where the Pritzsch ansatz is labelled as pattern (A). It is 
apparent that these six patterns are structurally different from one another. The question 
is whether their predictions for neutrino masses, flavor mixing angles and CP violation are 
distinguishable or not. 

The purpose of this paper is to answer the above question and to confront those six-zero 
textures of lepton mass matrices with the latest experimental data. Pirst, we shall present 
a concise analysis of the lepton mass matrices in Table 1 and reveal their isomeric features 
- namely, they have the same phenomenological consequences, although their structures are 
apparently different. Second, we shall examine the predictions of these lepton mass matrices 
by comparing them with the 2(T and 3cr intervals of two neutrino mass-squared differences 
and three lepton flavor mixing angles [10], which are obtained from a global analysis of 
the latest solar, atmospheric, reactor (KamLAND and CHOOZ [11]) and accelerator (K2K) 
neutrino data. We flnd no parameter space allowed for six isomeric lepton mass matrices 
at the 2(7 level. At the Scr level, however, their results for neutrino masses and lepton 
flavor mixing angles can be compatible with current data. Third, we incorporate the seesaw 
mechanism and the Fukugita-Tanimoto-Yanagida hypothesis [9] in the charged lepton and 
Dirac neutrino mass matrices with six texture zeros. It turns out that their predictions, 
including $23 ~ 45°, are in good agreement with the present experimental data even at the 
2(7 level. 



Because M/ and Mj/ are taken to be symmetric, a pair of off-diagonal texture zeros in M; or M( 
have been counted as one zero. 
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Let us begin with the diagonahzation of M; and Mi, hsted in Table 1. Without loss 
of generality, one may take their diagonal non-vanishing elements to be real and positive. 
Then only the off-diagonal non-vanishing elements of M/ and Mj, are complex. Each mass 
matrix M consists of two phase parameters (0 and (p) and three real and positive parameters 
(^4, B and C), as shown in Table 1, where their subscript "f or 'V" has been omitted for 
simplicity. The diagonahzation of M requires the following unitary transformation, 

/Ai 
U^MU* = As 

V A3 

where Aj (for i = 1,2, 3) denote the physical masses of charged leptons (i.e., Ai,2,3 
neutrinos (i.e., A^ = rrii). Due to the particular texture of M, U can be written as a product 
of a diagonal phase matrix (dependent on and (p) and a unitary matrix (independent of 
and ip), as illustrated by Table 1. The real parameters (^4, B, C) in M and (a^, h^, in U 
are simple functions of Aji 
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where a; = A1/A2 and y = A2/A3 have been defined. Note that a^, ^2 ^2 imaginary, 
and their nontrivial phases arise from a minus sign of the determinant of M(i.e., Det(M) = 
— 74C^e^*^). Since the charged lepton masses have precisely been measured [12], we have 
Xi 0.00484 and Di ^ 0.0594. On the other hand, < a;,^ < 1 is required by the solar 
neutrino oscillation data [1]. Hence < y,^ < 1 must hold, in agreement with Eq. (4). 
This observation implies that the isomeric lepton mass matrices under discussion guarantee 
a normal neutrino mass spectrum. 

The lepton flavor mixing matrix V, which links the neutrino mass eigenstates {vi, 1^2, 1^3) 
to the neutrino flavor eigenstates (z/g, t/^, 1/,-), results from the mismatch between the diago- 
nalization of M; and that of M^,. Taking account of Eq. (2), we obtain V — UlUi,, whose 
nine matrix elements read explicitly as 

Vv, = + (6y*6^e^^ + (c^^c^ , (5) 

where the subscripts p and q run respectively over (e,/x, r) and (1,2,3), and the phase 
parameters a and j3 are defined by a = {(p^^ — fi) ^ P and f3 = {(f) ^ ~ 'Pi)- It is worth 
remarking that Eq. (5) is universally valid for all six patterns of lepton mass matrices in 
Table 1. Hence they must have the same phenomenological consequences and can be referred 
to as the isomeric lepton mass matrices. 

Obviously, V consists of four unknown parameters: x^,, y^,, a and Their magnitudes 
can be constrained by current experimental data on neutrino oscillations. For the sake of 
convenience, we adopt the standard parametrization of V [13]: 



C12C13 S12C13 

V = -C12S23S13 - Sl2C23e"'^ -S12S23S13 + Ci2C23e"*^ 
V -C12C23S13 + Sl2S23e"*'^ -Sl2C23'Sl3 - Ci2S23e"*'^ 
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where = cos% and = sin^y (for ij = 12, 23, 13). Table 2 is a summary of the allowed 



ranges of two neutrino mass-squared differences (Amgi = 1712 — ml and Am|i = ml — m? 



31 — "^3 "n 



and three flavor mixing angles (sin^ 612, sin^ ^^23 and sin^ ^^13), obtained from a gobal analysis 
of the latest solar, atmospheric, reactor and accelerator neutrino data [10]. Because 



Am2i 2 ^ — X? 



R.^-^ = yt z f-2 (7) 

Am|i 1 - xiyl 



and 



sin^ 


^12 — 


\Ve2\' 


l-\Ve3\' 


sin^ 


6*23 = 




1 - IK3I' 


sin^ 


^13 = 


1^31^ 



(8) 

are all dependent on x,y, y^, a and /?, the latter can then be constrained by using the 
experimental data in Table 2. Once the parameter space of {x^^y^,) and (a,/?) is flxed, one 
may quantitatively determine the CP- violating phases (5, p, a) and the Jarlskog invariant J 
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(= Im[T42V^3K3^/r2]) for example [14]), which measures the strength of CP and T violation 
in neutrino oscillations. It is also possible to determine the neutrino mass spectrum and 
the effective masses of the tritium beta decay ((m)e = milKip + ^2|K2p + '^alKsP) and 
the neutrinoless double beta decay ((m)ee = l^iK^i + ^^2^2 + ^sKlD- ^he results of our 
numerical calculations are summarized as follows. 

(1) We find that the parameter space of {x,y,y„) or will be empty, if the best-fit 
values or the 2a intervals of Aml^, Amg^, sin^ 612, sin^ 623 and sin^ ^13 are taken into account. 
This situation is caused by the conflict between the largeness of sin^ 6*23 and the smallness 
of R^, which cannot simultaneously be achieved from Mi and M^, at the 2a level. 

(2) If the 3a intervals of Am|j, Amg^, sin^ ^12, sin^ ^23 and sin^ ^13 are taken into account, 
however, the consequences of Mi and M,y on neutrino masses and flavor mixing angles can 
be compatible with current experimental data. Fig. 1 shows the allowed parameter space of 
{xu,yu) and at the 3a level. We see that /5 ~ tt holds. This result is consistent with 
the previous observation [8]. Because of ~ 0.25, 777.3 ~ is a good approximation. 
The neutrino mass spectrum can actually be determined to an acceptable degree of accuracy: 
7773 ~ (3.8 - 6.1) X 10-2 eV, 7772 ^ (0.95 - 1.5) x 10"^ eV and 7771 ^ (2.6 - 3.4) x 10-^ eV, 
where x„ ^ 1/3 and yu ^ 1/4 have typically be taken. A straightforward calculation yields 
(m)e ~ 10^2 eV for the tritium beta decay and (m)ee ~ 10^^ eV for the neutrinoless double 
beta decay. Both of them are too small to be experimentally accessible in the foreseeable 
future. 

(3) Fig. 2 shows the outputs of sin^ ^12, sin^ ^23 and sin^ ^13 versus Ri, at the 3a level. It is 
obvious that the maximal atmospheric neutrino mixing (i.e., sin2^23 ~ 0.5 or sin2 2^23 ~ 1) 
cannot be achieved from the isomeric lepton mass matrices under consideration. We see 
that sin2 6'23 < 0.40 (or sin^ 26*23 < 0.96) holds in our ansatz, and it is impossible to get a 
larger value of sin^ ^23 even if Ri, approaches its upper bound. In contrast, the output of 
sin^ 9-12 is favorable and has less dependence on R,^. One can also see that only small values 
of sin2^i3 (< 0.016) are favored. More precise data on sin2^23, sin^ 6*13 and Ri, will allow 
us to check whether those isomeric lepton mass matrices with six texture zeros can really 
survive the experimental test or not. 

(4) We calculate the CP- violating phases (5, p, a) and the Jarlskog invariant J, and 
illustrate their results in Fig. 3. The maximal magnitude of J is close to 0.015 around 
S ~ 37r/4 or 57r/4. As for the Majorana phases p and a, the relation {p — a) ~ 7r/2 
holds. This result is attributed to the fact that the matrix elements (03, &2) '^2) '^^ are 
all imaginary and they give rise to an irremovable phase shift between Vpi and Vp2 (for 
p — e,iJ,, t) elements through Eq. (5). Such a phase difference may affect the effective mass 
of the neutrinoless double beta decay, but it has nothing to do with CP violation in neutrino 
oscillations. 

We proceed to discuss a simple way to avoid the potential tension between the smallness 
of Ri, and the largeness of sin^ ^23 arising from the above isomeric lepton mass matrices. In 
this connection, we take account of the Fukugita-Tanimoto-Yanagida hypothesis [9] together 
with the seesaw mechanism [15] - namely, the charged lepton mass matrix M; and the Dirac 
neutrino mass matrix Mq may take one of the six patterns illustrated in Table 1, while the 
right-handed Majorana neutrino mass matrix Mr takes the form Mr = MqI with Mq being 
a very large mass scale and I denoting the unity matrix. Then the effective (left-handed) 
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neutrino mass matrix Mi, reads as 



Ml 
Ma 



(9) 



For simplicity, we further assume Md to be real (i.e., = = 0). It turns out that the 
real orthogonal transformation Ub, which is defined to diagonalize Md, can simultaneously 
diagonalize M^,: 



U^MMd 
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mi \ 
m2 
may 



(10) 



where m^ = (if/Mo with di standing for the eigenvalues of Md. In terms of the neutrino 
mass ratios = mi/ma = {di/d'^f' and = 1712/1713 — (^2/(^3)^, we obtain the explicit 
expressions of nine matrix elements of Ui, — Ub' 
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(11) 



The lepton flavor mixing matrix V — UjUi, remains to take the same form as Eq. (5), but 
the relevant phase parameters are now defined as a = —ipi — (5 and P = — 0^ . Comparing 
between Eqs. (4) and (11), we immediately see that the magnitudes of (^^12, ^23) ^13) in the 
non-seesaw case can be reproduced in the seesaw case with much smaller values of x,^ and 
y,^. The latter will allow Ri, to be more strongly suppressed. It is therefore possible to 
relax the tension between the smallness of R^, and the largeness of sin ^23 appearing in 
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the non-seesaw case. A careful numerical analysis of six seesaw-modified patterns of the 
isomeric lepton mass matrices does support this observation. We summarize the results of 
our calculations as follows. 

(a) We find that the new ansatz are compatible very well with current neutrino oscillation 
data, even if the 2(t intervals of Am^^, Amg^, sin^^i2, sin^^23 and sin^ ^13 are taken into 
account. Hence it is unnecessary to do a similar analysis at the 3(7 level. The parameter 
space of {xi,,yu) and {a,P) is illustrated in Fig. 4, where Xiy ^ yu ^ 0.2 and /9 ~ vr hold 



approximately. Again ?a y Amli is a good approximation. The values of three neutrino 
masses read explicitly as ms (4.2 — 5.8) x 10"^ eV, m2 ~ (0.84 — 1.2) x 10~^ eV and 
mi fti (1.6 — 1.9) X 10~^ eV, which are obtained by taking Xy 0.2. It is easy to arrive 

at (m)e ~ 10~^ eV for the tritium beta decay and (m)ee ~ 10~^ eV for the neutrinoless 
double beta decay, thus both of them are too small to be experimentally accessible in the 
near future. 

(b) The outputs of sin^^i2, sin^ ^23 and sin^ ^13 versus are shown in Fig. 5 at the 
2a level. One can see that the magnitude of sin^ 6*12 is essentially unconstrained. Now the 
maximal atmospheric neutrino mixing (i.e., sin^^23 ~ 0.5 or sin^2^23 ~ 1) is achievable in 
the region of Ri, ~ 0.036 — 0.047. It is also possible to obtain sin^ ^13 < 0.035, just below 
the experimental upper bound [11]. If sin^2^^i3 > 0.02 really holds, the measurement of ^13 
should be realizable in a future reactor neutrino oscillation experiment [16]. 

(c) Fig. 6 illustrates the numerical results of (5, p, a and J. We see that \J\ ~ 0.025 
can be obtained. Such a size of CP violation is expected to be measured in the future long- 
baseline neutrino oscillation experiments. As for the Majorana phases p and cr, the relation 
a p holds. This result is easily understandable, because is real in the seesaw case. It is 
worth mentioning that the effective neutrino mass matrix does not persist in the simple 
texture as Mi has, thus the allowed ranges of 5, p and a become smaller in the seesaw case 
than in the non-seesaw case. 

Note that the eigenvalues of Md and the heavy Majorana mass scale Mq are not specified 
in the above analysis. But one may obtain \di/d2\ = ^/x^ ~ 0.4 and 1^2/^3! = ^Jy^ ~ 0.4. 
Such a weak hierarchy of (|c?i|, \d2\^ \d'A\) means that Md cannot directly be connected to the 
charged lepton mass matrix M;, nor can it be related to the up- type quark mass matrix (Mu) 
or its down-type counterpart (Mj) in a simple way. If the hypothesis Mr = MqI is rejected 
but the result U'^MyUy — Diagjmi, m2, ^3} with given by Eq. (11) is maintained, it 
will be possible to determine the pattern of Mr by means of the inverted seesaw formula 
Mr = Mj^Mj7^Md [17] and by assuming a specific relation between M^ and M^. For 
example, one may simply assume Md = with M^ taking the approximate Fritzsch form. 



Just for the purpose of illustration, we typically input Xp ^ y^ ^ 0.18 as well as mu/rric ~ 
rric/mt ~ 0.0031 and nit ~ 175 GeV at the electroweak scale [18]. Then we arrive at 
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in unit of GeV. This order-of-magnitudc estimate shows that the scale of Mr is close to that 
of grand unified theories Aqut ~ 10^^ GeV, but the texture of Mr and that of Mq (or Mi) 
have little similarity. It is certainly a very nontrivial task to combine the seesaw mechanism 
and those phenomenologically-favored patterns of lepton mass matrices. In this sense, the 
simple scenarios discussed in Ref . [9] and in the present paper may serve as a helpful example 
to give readers a ball-park feeling of the problem itself and possible solutions to it. 

In summary, we have analyzed six parallel patterns of lepton mass matrices with six 
texture zeros and demonstrated that their phenomenological consequences are exactly the 
same. Confronting the predictions of these isomeric lepton mass matrices with current 
neutrino oscillation data, we find that there is no parameter space at the 2a level. They 
can be compatible with the experimental data at the 3a level, but it is impossible to obtain 
the maximal atmospheric neutrino mixing. We have also discussed a very simple way to 
incorporate the seesaw mechanism in the charged lepton and Dirac neutrino mass matrices 
with six texture zeros. It is found that there is no problem to fit current data even at the 
2(7 level in the seesaw case. In particular, the maximal atmospheric neutrino mixing can 
naturally be reconciled with a relatively strong neutrino mass hierarchy. The results for 
the effective masses of the tritium beta decay and the neutrinoless double beta decay are 
too small to be experimentally accessible in both the seesaw and non-seesaw cases, but the 
strength of CP violation can reach the percent level and may be detectable in the future 
long-baseline neutrino oscillation experiments. 

We conclude that the peculiar feature of isomeric lepton mass matrices is very suggestive 
for model building. We therefore look forward to seeing whether such simple phenomeno- 
logical ansatze can survive the more stringent experimental test or not. 

One of us (S.Z.) is grateful to the theory division of IHEP for financial support and 
hospitality in Beijing. This work was supported in part by the National Natural Science 
Foundation of China. 
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TABLES 



TABLE I. The isomeric lepton mass matrices {Mi and Mj,) with six texture zeros and the 
unitary matrices {Ui and U,y) used to diagonahze them, where the subscripts "Z" and 'V have 
been omitted for simphcity. 
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TABLE IL The best-fit values, 2a and 3a intervals of Am^i, Am^i, sin^0i2, sin^ ^?23 and 
sin^ ^13 obtained from a global analysis of the latest solar, atmospheric, reactor and accelerator 
neutrino oscillation data [10]. 
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FIG. 1. The parameter space of {xi,,y^,) and (a,/3) at the 3(7 level. 
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FIG. 2. The outputs of sin^ 9i2, sin^ 623 and sin^ ^13 versus Ri, at the 3cr level. 
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FIG. 4. The parameter space of {Xu^y^j) and (a,/?) at the 2a level in the seesaw case. 
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FIG. 5. The outputs of sin^ $12, sin^ $23 and sin^ ^13 versus Rj, at the 2a level in the seesaw case. 
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FIG. 6. The outputs of {6, J) and (p, a) at the 2(T level in the seesaw case. 
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